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EXAMPLE OF TWO DIFFERENT POTENTIALS WHICH HAVE 
PRACTICALLY THE SAME FIXED-ENERGY PHASE SHIFTS 
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The Newton-Sabatier procedure for finding the potential from fixed-energy phase shifts is analyzed. A method is 
proposed for finding two quite different spherically-symmetric real-valued, piecewise-constant, compactly supported 
potentials which generate at a fixed energy the phase shifts which are practically indistinguishable for all £. In 
particular, an explicit concrete example of two such potentials qj, j = 1,2, is demonstrated. These potentials have 
the properties: 

1) sup |gi — 92! > 3, and qj, j = 1, 2, are of order of magnitude 1, 

2) 5p> = for i = 0, ... ,4 and < 10-5, Z > 4. 
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1 Introduction 



Finding real-valued potentials q{r), r = \x\, from the set {(5;};=o.i,... of the fixed-energy phase 
shifts (say at fc = 1), is an inverse problem which was much discussed in the literature 
[ 5 1 , 1^ , 1^] . Physicists worked much on this problem. They used a parameter fitting procedure 
[2| which consists of the following steps: 

Step 1. Given phase shifts Si at a fixed energy, which we take to be fc^ = 1 without loss of 
generality, one solves an infinite system of linear algebraic equations for some coefiicients c; 
(system (12.2.7) in ||^; the references to formulas with such numeration here and below are 
given alawys to the formulas from [2]). 

Step 2. Given these coefficients c;, one constructs some symmetric kernel ((12.2.1) in and solves 
linear integral equation (12.1.2) for the kernel K{r, p). If K{r^ p) is found, then the potential 
q{r) is calculated by the formula (12.2.11): 

-?^^L(!1I)^ (1.1, 

r ar r 

This procedure is presented in [5] and [2] as an inversion method and many physicists consider 

it as one. In this introduction we analyze this procedure. 

Let us define an inversion method as a method which satisfies the following conditions: 

1) given the exact data, corresponding to a q{r) G Q, where Q is a class of potentials in 

which the inverse scattering problem is uniquely solvable, one can carry through theoretically (and 

not necessarily computationally) each step of the method (we neglect here the numerical errors 

unavoidable in the numerical implementation of any method), 
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2) the method yields the unique solution q G Q to the inverse problem. 

In this paper we take the class Q of compactly supported real-valued potentials q{r) which 

satisfy the conditions: \q{r)\'^r'^dr < oo, q{r) = for r > a, where a > is an arbitrary large 
fixed number which is called the range of the potential. For q G Q the inverse scattering problem 
with fixed-energy phase shifts as the scattering data has unique solution [7]. Let us assume that 
the given scattering data come from a potential in this class. The question is: can one reconstruct 
such a potential using the Newton-Sabatier (NS) procedure and assuming that this procedure can 
be carried through? We prove that the answer to the above question is no. 

Solving the inverse scattering problem includes the requirement that the potential obtained 
by an inversion method produces the original data from which it was reconstructed. In [2] this 
property is called consistency of the method. 

The basic results of our paper arc: 

1) it is shown that if one applies the NS procedure to the data generated by a potential q G Q, 
then, assuming that this procedure can be carried through, it yields the potentials q{r) which are 
analytic in a neighborhood of the positive semi-axis of the complex r— plane and may have a simple 
pole at the origin, and therefore is not equal to the potential q € Q with which we started. 

2) it is shown that the stability estimate ( see estimate (1.5) below) from [8] is practically accu- 
rate: an example is constructed of two different piecewise-constant compactly supported potentials 
which produce practically the same (up to the error of order 10~^ ) fixed-energy phase shifts for 

all e. 

Therefore even if both steps in the NS method with exact data, corresponding to a g S Q, 
can be carried through in the above class Q, they lead to a potential qi{r) ^ Q, which has very 
special property (P), namely the function rq{r) has to be analytic in a neighborhood of the positive 
semiaxis of the complex r— plane. This property the potentials from class Q do not have. 

Consistency of the NS procedure is not proved in [2] and [5] : assuming that the NS procedure 
can be carried through for the given data, there is no guarantee that the potential generated by 
this procedure produces the same data with which one started. 

In [2, p. 205] an argument is given to prove the consistency of the NS procedure, but this 
argument is not completely convincing because: 1) it is not quite clear why A'^ and S'^ solve equation 
(12.2.5), 2) existence and imiqueness of the solution to Goursat-type problem (11.2.5)-(11.2.7) 
are not proved in [2]; since the coefficients in this problem are variable and have non-integrable 
singularity at the origin, the desired uniqueness and existence do not follow from the known results; 
in the literature one can find results on the existence and uniqueness of the transformation kernels 
for equation (11.1.2), but these kernels depend on i, while the transformation kernel used in [2] 
and [5] should be independent on £. 

Note that the basic ansatz of the NS procedure says that the transformation kernel (12.2.3) is 
generated by the symmetric kernel (12.2.1) and condition (12.2.2) holds. 

There is no guarantee that this ansatz can yield (by solving system (12.2.4) (^nd then using 
formula (12.2.11)) a potential which can produce generic scattering data. Indeed, as we prove 
below, the potential which is obtained by the NS procedure (in the case when this procedure can be 
carried through) has to satisfy property (P), which generic potentials do not have. Therefore the 
usage of this ansatz is not fully justified: there is no guarantee that one can fit the exact phase 
shifts, generated by a potential q G Q, using the potentials produced by the NS procedure based on 
the above ansatz. 

In the discussion of these principal points we are not talking at all about the difficulties which 
come from the fact that in practice the data are never exact. 

There is no guarantee that the steps 1 and 2 of the NS procedure with exact data can be carried 
through, in general, because: 
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1) System (12.2.7) may have, in general, no solution in the class (12.2.2), 

2) Integral equation (12.1.2) for the kernel K{r, p) may be not solvable for some (not too small) 
r > 0. 

Note that existence of the transformation kernel (11.2.4) is not proved in [2] and [5] but was 
proved recently in [11]. 

It is easy to prove that equation (12.1.2) is imiquely solvable for sufficiently small r > since 
in this case the norm of the corresponding integral operator in C(0, r) tends to as r ^ 0. 

If it is possible to solve (12.2.7) in the class (12.2.2) and equation (12.1.2) is uniquely solvable 
for all r > 0, then the potential (12.2.11) must have the following special property: 

(P) The function rq{r) must be an analytic function of r in a neighborhood of the positive 
semiaxis (0, oo). 

The proof of property (P) for the potentials obtained by the NS scheme is given below. 
From the logical point of view property (P) is sufficient for the proof of our basic conclusion: 
Conclusion: For the generic data corresponding to a compactly supported potential q{r) G Q 
the NS procedure cannot reproduce the original potential. 
Thus, we have the following logical possibilities: 

i) Either the NS procedure can be carried through for the phase shifts corresponding to a q G Q 

and the reconstructed by this procedure potential qi{r) reproduces the original phase shifts, and 
then the inverse scattering problem has more than one solution and the solution q € Q cannot be 
obtained by the NS procedure; 

ii) or the NS procedure cannot be carried through for the phase shifts corresponding to q G Q; 
Hi) or the NS procedure can be carried through for these phase shifts but the obtained potential 

qi (r) does not generate the original phase shifts, in which case the NS procedure does not solve the 
inverse scattering problem. 

Similar argument holds for any class of potentials which do not have property (P) and for which 
the uniqueness of the solution to the inverse scattering problem with exact data is valid. 

Let us now prove that the potential, obtained by the NS procedure, assuming that this procedure 
can be carried through for the given data, must have property (P). 

Proof of property (P). 

The bmctions Uf in (12.2.3) arc entire functions of r which have the following known asymptotics 
for large £: U£ ~ \/1"(2||t)~^(2^ + 1)^^/^. Therefore the function (12.2.1), with q satisfying 
(12.2.2), is an entire function of r and r' . 

Let us first prove that rq{r) is analytic in a neighborhood of the origin, since this claim already 
is sufficient for a justification of our basic conclusion and the proof of this claim is quite easy. 

The system (12.2.4) has the matrix i«'(r) = ue{p)ue'{p)p~'^dp (see formula (11.4.5) with 
Vq = 0). Since Ui{r) is an entire function of r and is 0{r^^^) as r ^ 0, one concludes that Lu'{r) 
is an entire function of r. For sufficiently small r system (12.2.4) has dominating diagonal terms 
and therefore is uniquely solvable by iterations. Since U(,{r) and Lu'{r) are entire functions of 
r, the Neumann scries yields the solution ^i{r) of (12.2.4) which is an analytic function of r in 
a neighborhood of the point r = 0. This implies that the kernel K{r,r'), defined by (12.2.3), is 
analytic with respect to both variables in a neighborhood of r = 0, r' = 0. Therefore K{r, r) is 
analytic in a neighborhood of r = and consequently the potential (12.2.11) (defined by formula 
(1.1)) is analytic in a punctured at r = neighborhood of the origin and the function rq{r) is 
analytic in a neighborhood of the origin. Our first claim is proved. 

Let us now prove that rq{r) is analytic in a neighborhood of the positive semiaxis provided 
that equation (12.1.2) is uniquely solvable. The kernel and the free term of this equation, as we 
have proved above, are entire functions of r and r' . It is known that the solution of a Predholm 
second kind integral equation whose kernel and free term depend on a parameter z analytically in 
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a certain region A of a complex plane, is analytic in a neighborhood of a point zo £ A, provided 
that the above Fredholm equation is uniquely solvable at z — zq (see [10] for a more general result). 

Let us introduce the new variables p — rs,r' = rt in equation (12.1.2) with the aim to get 
equation with a fixed integration region. Equation (12.1.2) becomes: 
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K{r,rt) = f(r,rt) — j K{r,rs)f{rs,rt)r s ds. 

Let 







a{s,t,r) := f{rs,rt){rst) \ K{r,rs)s ^:=6(s,r), f{r,rt)t ^ -.^ g{r,t). 
Then the above equation takes the form: 

b{t,r) = g{r,t) — / a{s,t,r)b{s,r)ds. 
Jo 

This equation is equivalent to equation (12.1.2) in [2]. It is a Fredholm second kind equation 
whose kernel a{s,t,r) and the free term g(r,t) depend on the parameter r analytically on the 
whole complex r— plane because the kernel f{s,t) is an entire function of s and t, as we proved 
above, and it vanishes at s = and at t = 0. Also g{r, t) and a(s, t, r) are entire functions of t. If 
the NS procedure can be carried through, then this equation has to be uniquely solvable for any 
r > 0. Therefore its solution is an analytic function of r in a neighborhood of positive semiaxis 
(0, oo) and an analytic function of t on the whole complex plane. This implies that the kernel 
K{r,r) is an analytic function of r in a neighborhood of the positive semiaxis (0, oo). Thus, the 
potential (12.2.11) has property (P). The proof is complete. □ 

Therefore, in general, the inverse scattering problem is not solvable in the class Q by means of 
the NS procedure. 

The uniqueness of the recovery of q{r) from the fixed-energy phase shifts (5;, Z > 0, is not 
established in ||]. Such a uniqueness theorem follows for g G Q from the uniqueness theorem 
proved in for potentials which are not necessarily spherically symmetric. 

No such theorem is known for potentials which decay as some power of -, that is. 



\l\ = 0[j^^), m>2isfixed. (1.2) 

It is known (see [|j) that the range a of a compactly supported potential, which does not change 
sign in some neighborhood (a — s,a), can be calculated in terms of {Si} : 

a= hm C^tLllSilii] . (1.3) 
i-^oc \ e / 



It follows from (1.3) that 



f ae\l + o(l)]Y^ , , 

-''I 27Tl ) ' '^ = ^°"'^^>0. (1.4) 

This implies that (5; decay very fast as I — > cx), provided that q(r) is compactly supported. 

If q(r) decays on a power scale (L2) then bi decays also as some positive power of (see [||). 

In this paper we construct two piecewise-constant, real- valued, compactly supported potentials 
gi(r) and <?2(''), which are quite different (see Table 1 in sec. 3) and which generate practically 
identical phase shifts Si. 
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The potentials qj, j = 1,2, are constructed by a computer code, which recovers q{r) (in the 
above class of potentials) from the data {Si}, I > 0. 

It is known that the inverse problem of quantum scattering with fixed-energy scattering data 
is very ill-posed, if the energy is not very high, and the numerical reconstruction of the potential 
from noisy scattering data is quite difficult . 

However, our numerical results also show that if an initial approximation to the unknown 
potential is known with relatively small error (« 10""'^), then the recovery of q(r) is possible with 
high accuracy (10^^''). 

In ^ a mathematically justified method is proposed for stable inversion of fixed-energy scat- 
tering data for not necessarily spherically symmetric potentials with compact support. If the 
potential q{r) is bounded spherically , symmetric and has compact support then the result in ||^ 
yields an algorithm for reconstruction of q(r) from the phase shifts Si known with some error e. 
This algorithm yields qe such that 

, . ^ , (ln|lne|)2 ^ ^ 

supm -qe <c ' , e^O, 1.5 

? |lne| 

where c > is some constant which does not depend on e, but does depend on some norm of 



the potential, (j(^) is the Fourier transform of q{x). The right-hand side of (L5) tends to zero 
as £ — > 0, however it tends very slowly. Therefore practically one can hope to recover q{r) from 
the fixed-energy scattering data, that is, from the phase shifts 6i at a fixed k, only if either the 
class of the potentials is a priori chosen to be rather narrow, for instance, finite-parametric, or the 
accuracy of the data is very high. Our numerical results confirm these conclusions. The example 
we constructed is an illustration of the stability estimate (1.5). 

In section 2 a numerical method is described which we use to invert the fixed-energy phase 
shifts for the potential q{r). In section 3 the numerical results are presented. 



2 Numerical method 



Consider a finite set of points — < ri < r2 < ... < — R and a piecewise-constant 
potential 

lij") = on [^i-iifi) for i = 1, . . . , TV, and g = for r > R. (2.1) 

Denote k| := k^ — qi, where i = 1,...,N, and k is some fixed positive number. Consider the 
following problem for the radial Schrodinger equation: 

^+{k^-'-^^y = m^ lim(2; + l)!!r--V.W = l, (2-2) 
which we rewrite as: 

^+(-rV^j^.=0 (2.3) 

on the interval r^^i < r < r^. On [ri^i,ri) one has the following general solution of ( |2.3| ) 

ipi{r) = AiUi{Kir) + BiVi{Kir), (2.4) 



where 



Ul{z) = \l'^Jl+l/2[z), Vi{z) = J'^Yl+l/2{z). 
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We assume below that Ki does not vanish for all i. li Hi — Q for some i then our approach is still 
valid with obvious changes. 

From the regularity of ipi at zero one gets Bi = 0. Denote Xi — Bi/Ai, then xi = 0. We 
are looking for the continuously differentiable solution ipi. Thus, the following interface conditions 
hold: 

AiUi{Kiri) + BiVii^KiTi) = Ai+iUi{Ki+iri) + Bi+iVi{Ki+iri), 



■jT^lAiu'iiKin) + Biv'iiKiri)] = Ai+iu'i{Ki+iri) + _B^+l^;,'(K^+l^i). 
The Wronskian W{ui{r),vi{r)) = 1, thus 

Ai+i = v'i{Ki+iri)[AiUi{Kin) + BiVi{Kiri)] - -^vi{Ki+in)[Aiu'i{Kiri) + BivKK^ri)], 

Bt.+i = -;T^uiiKi+in)[Aiu'i{K^ri) + Biv'i{Kiri)] ~ u'i{Ki+iri)[AiUi{Kiri) + BiVi{Kin)]. 
Therefore 

^i+l \ _ ^ f "ll "12 \ f 



where the entries of the matrix a' can be written explicitly: 







- KiU[{Kin)vi{Ki+iri 


"12 






"21 




Ki+lUl{Kiri)u'i{Kt+m 


"22 







(2.5) 



(2.6) 



(2.7) 



(2.; 



Thus 

x.+i = i^ + "P% xr.= ^ (2.9) 

"ll+"l22;i 

Denote the phase shift S{1, k) by the formula 

^-^^ sin(fcr - ^ + 6{l, k)) r ^ 00 , (2.10) 
where F{1, k) is the Jost function. For r > R 

(pi{r) = AN+iUi{kr) + BN+iVi{kr). (2.11) 



From (2.11) and the asymptotics uiikr) ^ sin(fcr — 17:12), vi{kr) ~ — cos(/cr — It:/ 2), r 00, one 
gets: 

ta,iid{k,l) — — = —xn+1, 6(k, I) = — arctanxN+i- (2-12) 

An+1 

Now let us assume that the column vector S* = (Sq, . . . , 6lf_i) is given and introduce the nonlinear 
operator $ : — > R^ , which maps a column vector k = (ki, . . . , kn) into a vector $(k) — 5 — S* , 
where 6 is given by ( 2.12 ) and S* , j — 0, . . . , N — 1, are given fixed energy phase shifts. Thus, 
we are looking for the solution to the problem 

$(k) = 0. (2.13) 
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In order to solve this equation by means of the Newton method, first one has to find the Jacobian 
||$'|| = \\^j<^>^\\, where dj = d/dn^. 

Since a' — a'^^Hi, k^+i), Xi = Xi^ni, . . . , Ki), one has 



deta* 
("ii + "122; 

and, for j = i or j = i + l, one gets: 



(ail + c(i2Xi){djali + a;»gj-Q;|2) " ("21 + Q^22'^»)('9j«ii + a^^^j-ala) 

("11 + ai22;i)^ 

where 

aiOfli = Ki+iriu'i{Kiri)v'i{Ki+iri) - u'i{Kiri)vi{ni+iri) - Kinu'/ {Kiri)viiKi+in), 

dia'22 = v'i{Kin)ui{Ki+in) + KiTiV'/ {Kiri)ui{Ki+in) - Ki+inv'i{Kiri)u[{Ki+iri), 







di+ia{2 




di+ia'21 


— ^i^i'^i 


di+ia'22 





Then one solves problem ( 2.f 3 ) iteratively: 



J = 0,1,..., ■.= {^{,....,n'^), (2.16) 



(2.15) 



where the step 7 is chosen small enough. One stops the iterative process when ||$(k^)|| becomes 
smaller than some small threshold e. 



Instead of inverting matrix $' in ( 2.16 ) one denotes h-' := [<I>'(k^)]~^$(k-') and uses the Gauss 
process with the choice of the maximal element by lines and columns of the corresponding linear 
system: 

[<S?'{K.^)]i^^hl = K^+i = - jy. (2.17) 

If at a step Iq of the Gauss process, the maximal element of the remaining elements of the matrix 
becomes smaller than some threshold ei one sets all the remaining quantities h", Z > equal to 
zero. 

3 Numerical results 



7 



As was pointed out in the Introduction, the principal difficulty of the numerical realization of 
the inversion schemes, developed for the fixed energy inverse scattering problem, is the fast decay 
of the phase shifts Si . Therefore one cannot discriminate between the noise and the data starting 
from relatively small I. So one can use for the numerical reconstruction of a potential very few 
phase shifts Sq,. . .,Si„, thus neglecting the rest, which we call the "tail" Si, I > Iq. The goal of 
our numerical experiments is to show numerically that in general the " tail" of the data can not be 
neglected. 



In the numerical experiments we start with some piecewise-constant potential (2.1) and take 
N ~ 10. For example let us start with the following potential: 

q^ = {l + cos{i/2))e-\ i = 1, 2, . . . , 10. (3.1) 



This potential is denoted in Table 1 by qorig- Then, using formulas (U), (|2T|) we calculate the 



phase shifts Si, I = 0,1,..., 9, corresponding to some fixed value of k {k=2). Then taking the 
first five phase shifts Sq,.. .,S4 unchanged, we perturb the phase shifts S^,.. .,Sq a little, by a 
random quantity of of order of magnitude ~ 10~^. This corresponds to the order of ^5. By the 
numerical method described in the previous section (with ei = 10"^" and 7 = 1) we reconstruct 
the piecewise-constant potential (denoted in Table 1 as qrec), which generates the perturbed phase 
shifts Si with the discrepancy e = 10~^^. The corresponding numerical results are presented in 
Table 1 below. One can see that the first five phase shifts Si and Si of qorig and qrec correspondingly 
coincide, the remaining five phase shifts differ by the quality of order 10~^, but the potentials qorig 
and qrec are quite different. 



Table 1. 



n 


Qorig 


Si 


5i 


qrec 


0.5 


0.6907240 




-9.941752 


10 


-2 


-9.941752 


10 


-2 


-2.415259 


1 


0.2084572 




-3.779873 


10- 


-2 


-3.779873 


10 


-2 


1.558406 


1.5 


5.330886-10 


-2 


-1.179639 


10- 


-2 


-1.179639 


10 


-2 


-0.589802 


2 


1.069364-10 


-2 


-3.014222 


10 


-3 


-3.014222 


10 


-3 


0.355841 


2.5 


1.339883-10 


-3 


-6.494566 


10- 


-4 


-6.494566 


10- 


-4 


-0.171777 


3 


2.480612 - 10 


-5 


-1.691416 


10- 


-4 


-1.719357 


10- 


-4 


8.157301-10-2 


3.5 


5.794400 - 10 


-5 


-1.026825 


10- 


-4 


-9.611266 


10 


-5 


-3.368591-10-2 


4 


1.161896- 10" 


-4 


-7.547580 


10- 


-5 


-6.558222 


10 


-5 


1.191644-10-2 


4.5 


9.739553-10 


-5 


-4.157540 


10- 


-5 


-3.745421 


10 


-5 


-3.106548-10-3 


5 


5.827817- 10" 


-5 


-1.675351 


10 


-5 


-2.219373 


10- 


-5 


5.785574-10-'* 
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